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1 Introduction 



In the field of Inertial Confinement fusion(ICF), laser-produced plasmas generally fall be- 
tween collisionless and collisional regimes and could not simply be considered as being in 
local thermal equilibrium because of small size and short duration of the plasmas. 

When a plasma is weakly coupled, i.e. the Coulomb logarithm InA of the plasma is much 
larger than 1, transport process is mainly governed by the multiple small- angle Coulomb 
scattering between charged particles, and can be accurately described by Fokker-Planck(FP) 
equation. By comparison, a classical model named Spitzer-Harm theory has been proposed 
for long time. However, steep electron temperature gradient around critical density make 
the classical Spitzer-Harm theory failure to describe the electron thermal conduction in this 
region. Since energetic electrons carrying the heat flux have a long mean-free path that 
in the same order of or even longer than the scale length of the electron temperature, the 
electron heat flux is determined by the electron temperature distribution rather than the 
local temperature gradient. The electron transport becomes non-local. 

Much of the progress in the Fokker-Planck simulation of a plasma has been possible follow- 
ing the pioneer work of CHANG and COOPER. Based upon previous work, EPPERLEIN 
has advanced an implicit and conservative difference scheme for the Fokker-Planck equation. 
In this work, we apply Eppcrlcin's scheme to develop a Fokker-Planck code written in C 
language by Doc. Bin Zhao. 



2 Fokker-Planck equations and numerical scheme 

The electron Fokker-Planck equation for a high-Z laser produced plasma can be written as: 

% + v -ji--— E -% = c M + c ^(f) (i) 

at or m e dv 

C ee is the electron-electron collision operator, and C e i is the Lorentz electron-ion collision 
operator. Then distribution function f(x,v,t) can be expanded in a series of Legcndrc 
polynomials: 

oo 

f(x,v,t) = J2fl(x,v,t)P l (n) (2) 

1=0 

where yu = v x /v = cosO is the direction cosine. In the diffusive approximation, the series 
(2) is truncated at N = 1, namely, only /o and f\ are kept in the expansion, and the time 
derivative of /i is neglected. It is easy to obtain the equation for the isotropic part of the 
distribution function which is coupled with the anisotropic part f\\ 

dfo vdh _ 1 8 ( v 2 eE , 

For high-Z plasma, the anisotropic part of the electron-electron collision term is negligible 
in comparison with that of electron-ion collision term. The equation for fx can be simplified 
as: 



^ = 1 ,df eEdfp. 

v e i dx m e dv ' 



where v ei — 4:irN e Ze 4 /m 2 lnA/v 3 , Z is the ion charge state, N e is the electron number 
density. Then we obtain a closed equation of fo- 

dfo = d_r dfp a dfp -, Id 

dv dv L dx v dv * v 2 dv ^ 



[ x (a 2 ^)-av^]+C ee (fo,fo) + S IB 
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where a = eE/m e .x = v 2 /3f e i- The electron-electron collision operator C ee (/o,/o) is given 
by 

C ee (/o,/o) = T ee ^[C(f )fo + D(f )^} (6) 

where T ee — Aire 4 /m 2 lnK — v^/N e T,T is the electron-electron collision timc.iv is the electron 
thermal velocity, and C(/o) and D(f ) are defined as 



•> = /"' 

Jo 



C(/ ) - / duu 2 f (u,t) (7) 



i i-v 2 /-+00 



duuf (u,t) (8) 

The self-consistent eletrostatic field -E or equivalently a can be obtained from the quasi- 
neutrality condition. Taking the zeroth order velocity moment, we obtain the continuity 
equation for the electron density N e . 

dt dv^ 

where T is the electron flux. 



r = -4,fV<M%-?f) 

Jo OX V ov ' 



(9) 



The condition of quasi-neutrality requires a divergence-free partial flux T. With the assump- 
tion of the immobile ion background, the self-consistent eletrostatic field a can be obtained 
with the condition r = 0, 

= fdvv 2 x(df /dx) 
a Jdvv 2 ( X /v)(dfo/dv) ( ' 

With the condition of quasi-neutrality, the second term on the right hand side of Eq. (5), usually 
called Ohmic term, can be neglected. 

On the right hand side of Eq.(5), we add an additional term Sjb to include the IB absorption 
of the laser energy. According to LANGDON's theory, Sib can be written as 

*.(*)- ^^(if) 

In order to ensure better energy conservation, we calculate D directly from the equation 

— (vD)=v 2 J^ duuf (12) 

Phase-space differentiation is carried out with on a fixed Cartesian grid in (x, w).The spatial 
grid is usually uniform and the grid of the velocity space can be feathered with finer resolution 
at low velocity grid when the laser-heating effect is taken into consideration in the simulation. 
In velocity space, number conserving boundary conditions are adopted. In order to keep a 
zero heat flux in the spatial domain boundary, fixed boundary conditions are implemented 
in spatial space. 

My work mainly concentrates on the comparison between classical thermal transport by 
using Spitzer-Harm theory and the Fokker-Planck simulation. Some examples concerned 
with different aspects are shown in following sections. 
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3 When the initial density has linear profile 

3.1 Case 1: n e = 10 14 cm~ 3 

In this example of simulation, the conditions are as follows: 

• no laser heating. 

• density profile: linear from 4ev to 5ev, 4ev corresponds to the left point and 5ev 
corresponds to the right point. 

• only considering a 1-D model with the length L — 200cm, the space is divided into 
100 girds. 

• time step is St = 1.0 x 10~ 6 s, which means after St all the parameters of the system 
will be recalculated. After every 100 x St, the code generates an output. 

The simulation results are shown as follows (Figure l,t = 20000?is), the dark line represents 
the result by using Fokker-Planck equation, while the red one represents the result by using 
classical SH model. 




Figure 1: Heat Flux 



The mean free path is 0.256cto at this time point, which is much shorter than the width of 
each grid(2cm).The electron-electron collision time is 1.22 x 10~ 9 s. 

By examining the results, the two different models bring out nearly the same results. This 
indicates that under this certain conditions, the classical theory is still applicable. 
The following picture [Figure 2] plots the temperature profile at t = 20000?is. 

3.2 Case 2: n e = KFcm" 3 

In this example of the simulation, the mainly difference from case 1 is the electron density. 
Meanwhile, some simulation parameters have been changed: 
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Figure 2: Temperature Profile 



• the density is 10 11 cm~ 3 

• time step is St = 1.0~ 8 s 

The following picture [figure 3] shows the simulation results. From the picture, the two different 
results show distinct difference in calculating heat flux. This means, in this certain case, due 
to the low density of electron,the local model should be applied and the classical theory is 
not valid at all. 
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Figure 3: Heat Flux 
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The following picture [figure 4] shows the temperature distribution along x axis when t = 
2000ns. It is no doubt that the final temperature will be all uniform in space. 




Figure 4: Temperature Profile 



4 When the classical heat flux remains constant at ini- 
tial state 

When the classical heat flux remains constant, the initial temperature profile will be in a 
certain mathematical form. The classical thermal transport can be written as: 

Q sh = C ■ T 5 ' 2 VT (13) 

If heat flux remains a constant, say Co, then we can get: 

C ■ T 5 / 2 VT = C (14) 

Thus the temperature should be 

T = {C x x + C 2 ) 2/7 (15) 

So if given the temperature value of two points along the axis, the C\ and Ci can be solved. 

In the simulation, the number of grid is 400, and the number of the grids of velocity space is 
80. Here we consider a period as long as 400cm. We set the first 50 girds to be 4ev and the 
last 40 grids to be 5ev,in order to better analyze the evolution of temperature according to 
time. Like the above discussions, two examples of n e = 10 14 cm -3 and n e = 10 11 cm -3 . 

4.1 Case 1: n e = 10 14 cm~ 3 
As for this case, the parameters are listed below: 
• time step is St = 1.0 x 10 _8 s. 
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• for every 100 steps, the program generates an output. 
The result presented in Figure 5 when t — 2000ns. The result shows that it is a classical 
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Figure 5: Heat Flux 

case, in which classical Spitzcr-Harm theory can get the right physical result. The calculated 
mean free path for electrons is 0.256cm, and the collision time of electrons and electrons is 
1.22 x 10~ 9 . 

The following picture gives the temperature profile at this time point. 
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Figure 6: Temperature Profile 
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4.2 Case 2: 
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The mean free path in this case is 181.5cm. Such long mean free path leads to totally local 
thermal transport. The time step is set 1.0 x 10 -10 s and for every 200 steps, we get a result. 

Figure 7 gives the heat flux when t — 40ns. 
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Figure 7: Heat Flux 



As for the temperature profile, it is shown in Figure 8. 




Figure 8: Temperature Profile 
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5 Conclusion 



In summary, we develop a Fokker-Planck code to investigate the electron transport process 
in laser-produced plasma in relevance to inertial confinement fusion. The diffusive approx- 
imation is adopted to simplify the Fokker-Planck equation. The equations are numerically 
implemented by infinite difference method and solved implicity in time. By comparing the 
classical Spitzer-Harm theory and the Fokker-Planck equation, the physical picture of local 
transport and nonlocal transport becomes much clearer. When the temperature gradient 
is deep or the density of electrons is low, the effect of local thermal transport will become 
dominant which gives total different results. Always is the case that the nonlocal thermal 
heat flux is less than classical heat flux within the calculated region. However, by now we 
only talk about a simple example in which there is no laser heating, and this case is common 
in MCF when the magnetic field is so weak that it can be neglected. As for some complex 
conditions, further simulations and discussions will be the future work we will do. 
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